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Antiferromagnetic Ising Model on Inverse Perovskite Lattice 
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We study thermodynamic properties of an antiferromagnetic Ising model on the inverse per- 
ovskite lattice by using Monte Carlo simulations. The lattice structure is composed of corner- 
sharing octahedra and contains three-dimensional geometrical frustration in terms of magnetic 
interactions. The system with the nearest-neighbor interactions alone does not exhibit any 
phase transition, leading to a degenerate ground state with large residual entropy. The degen- 
eracy is lifted by an external magnetic field or by an anisotropy in the interactions. Depending 
on the anisotropy, they stabilize either a 3D ferrimagnetic state or a partially-disordered anti- 
ferromagnetic (PDAF) state with a dimensionality reduction to 2D. By the degeneracy- lifting 
perturbations, all the transition temperatures of these different ordered states continuously 
grow from zero, leaving an unusual zero-temperature critical point at the unperturbed point. 
Such a zero-temperature multicriticality is not observed in other frustrated structures such as 
face-centered cubic and pyrochlore. The transition to the PDAF state is represented by either 
the first- or second-order boundaries separated by tricritical lines, whereas the PDAF phase 
shows 1/3 magnetization plateaus. 

KEYWORDS: geometrical frustration, inverse perovskite, Ising model, residual entropy, magnetization 
plateau, partial disorder, dimensionality reduction, ferrimagnetism 



Geometrical frustration is one of the long-standing 
issues in condensed matter physics. The problem has 
been intensively investigated in antiferromagnetic (AF) 
spin systems. Frustration, in general, suppresses the for- 
mation of simple-minded long-range orders (LRO) and 
opens the possibility of interesting phenomena such as 
a degenerate ground-state manifold with an extensive 
residual entropy, complicated ordering structures, and 
spin-liquid states. 1 ) 

Even in one of the simplest models, namely the Ising 
model, frustration leads to intriguing properties. 1 )' 2 ) A 
prototypical example is the AF Ising model on a trian- 
gular lattice. The system does not show any LRO down 
to zero temperature, and the ground state has macro- 
scopic degeneracy. 3 ) 5 ) Among other frustrated lattices, 
3D models have also attracted considerable attention 
partly because certain real materials are available as 
their experimental counterparts. For instance, the face- 
centered cubic (fee) lattice has been studied for the 
problem of binary alloys and the pyrochlore lattice has 



been examined to describe magnetite 6 ) and spin-ice com- 
pounds. 7 ) 

In this paper, we discuss a new 3D frustrated Ising 
model, an inverse perovskite AF Ising model. The in- 
verse perovskite lattice consists of a simple cubic lat- 
tice of octahedra which share their corners, as shown in 
Fig. 1(a). The corner-sharing geometry reminds us of the 
well-studied pyrochlore lattice, which consists of corner- 
sharing tetrahedra. Alternatively, the lattice structure is 
viewed as an alternating stack of 2D square lattices with 
different lattice constants along the [001] direction (or 
the equivalent [100] or [010] direction), while it is a stack 
of Kagome lattices along the [111] direction (or the equiv- 
alent directions such as [111]). These stacks are similar 
to another well-studied lattice, the fee lattice, which is 
a stacking of 2D square (triangular) lattices along the 
[001] ([111]) direction. The inverse perovskite structure 
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Fig. 1. (a) Inverse perovskite lattice structure, (b) A unit cell 
with three sublattice sites. J and Jj_ are the interactions in the 
Hamiltonian in eq. (1). 



Fig. 2. Temperature dependences of (a) specific heat per site and 
(b) entropy per site at J = J± and h = 0. The dashed line in 
(b) denotes the value in the high-temperature limit, log 2. 
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is found in several materials such as the superconductor 
MgCNi3 8 ) and the negative thermal expansion material 
MnatCu^Ge^N. 9 ) 

In the present study, we investigate effects of frustra- 
tion on the thermodynamics of the AF Ising model on the 
inverse perovskite lattice. By employing a Monte Carlo 
(MC) method, we study phase transitions in this model 
by tuning the degree of frustration and obtain the phase 
diagram in a wide range of parameters. In addition to 
some similar behaviors to the fee and pyrochlore sys- 
tems, we find a peculiar critical point where all the tran- 
sition temperatures are suppressed down continuously to 
zero due to a characteristic structure of the degenerate 
ground-state manifold in the inverse perovskite system. 

Our model Hamiltonian is defined in the form 

n = jJ2 s i s J + J ±J2 SiSj-hJ2Si, (i) 

(ij)e\\ (i,j)e± i 

where Si = ±1 denotes the Ising variable and h is an ex- 
ternal magnetic field. In the inverse perovskite structure, 
three sites are present in a unit cell, as shown in Fig. 1(b), 
and Si in eq. (1) symbolically represents S u ^, where v 
denotes the unit cell and a = 1, 2, or 3 is the sublattice 
index. The summations of the interaction terms are taken 
over the nearest-neighbor sites on the inverse perovskite 
lattice in Fig. 1(a), and an anisotropy is introduced in 
the interactions J and J_l, as shown in Fig. 1(b): J is for 
pairs of S V: \ and S v > ^ spins within the same xy plane, 
and J± is for pairs including S v ^ spins. We consider both 
cases of J > J± > and J± > J > 0. 

We study thermodynamic properties of the model in 
eq. (1) by using MC simulations. To avoid the critical 
slowing down in the frustrated system, we employ the 
exchange MC method 10 ) by an exchange sampling taken 
after every 5 standard MC steps. Typically, we perform 
10 5 -10 7 exchanges in total for the system size L 3 (N s = 
3L 3 spins) up to L = 64 under the periodic boundary 
conditions. We set the Boltzmann constant &b = 1- 

First, we study the isotropic case J = J± at the mag- 
netic field h = 0. Figure 2(a) shows the temperature (T) 
dependence of the specific heat per spin plotted in the 
logarithmic scale of T. The specific heat is calculated 
by C(T) = ((H 2 ) - (H) 2 )/T 2 , where the bracket de- 
notes the thermal average. As shown in Fig. 2(a), the 
specific heat shows no singularity, except for a broad 
hump at T/J ~ 2 which is almost system-size indepen- 
dent. This indicates that the model does not exhibit any 
phase transition at finite T. The entropy is calculated 
by S(T) = E(T)/T - / 1/T E(T') d(l/T) + S(oo) , where 
E(T) = (H) is the internal energy and S(oo) = 7V s log2 
is the entropy in the high-T limit. The result is plotted 
in Fig. 2(b). With lowering T, the entropy gradually de- 
creases and approaches a finite value as T — > 0. The resid- 
ual entropy per spin is estimated as 0.2941(1) ~ 42% of 
log 2 (statistical errors in the last digit are shown in the 
parentheses, hereafter). This result means that macro- 
scopic degeneracy remains in the ground state. The ab- 
sence of a phase transition and the large residual entropy 
are found in other 3D frustrated AF Ising systems such 
as pyrochlore 7 ) and garnet, n ) and appears to be a com- 
mon feature of corner-sharing- type frustrations. 



When a magnetic field is applied, the macroscopic de- 
generacy is lifted and a phase transition takes place at 
finite T. Figure 3(a) shows T dependences of the spe- 
cific heat. The specific heat shows a sharp singularity 
for \h\ > 0, indicating a phase transition; for instance, 
a sharp peak at T/J ~ 1.1 (2.0) for h/J = 0.5 (4.0) 
as shown in Fig. 3(a). The small hump at T/J ~ 0.3 
for h/J = 0.5 comes from the disordered spins remain- 
ing in the low-T ordered state (see below). Below the 
transition temperature, the magnetization process shows 
plateaus at (m) — ±1/3 as shown in Fig. 3(b), where 
m = J2v iSv,i/N s . The plateau state shows a partially- 
disordered AF (PDAF) ordering, which consists of the 
stacking of AF-ordered square lattices and disordered 
spins sandwiched between AF layers, as shown schemat- 
ically in the left inset of Fig. 4. This dimensionality re- 
duction, i.e., 2D LRO in a 3D system, yields 2 L -fold de- 
generacy due to the twofold degeneracy of the AF or- 
dering pattern in each layer, leading to (3 x 2 L — 3)-fold 
degeneracy in total. The factor 3 comes from the stack- 
ing directions x, y, and z, and —3 subtracts the double 
counting. The schematic picture in Fig. 4 shows one of 
these degenerate states. The magnetic ordering and de- 
generacy are confirmed by calculating the 2D AF order 
parameter in each xy, yz, and zx layer. We summarize 
the critical temperature T c determined from these mag- 
netic properties in the phase diagram in Fig. 3(c). T c is 
finite at < \h\/J < 8, showing its maxima T c / J ~ 2 at 
h/J~ ±4. The transition is of first order, confirmed by 
a detailed analysis of the probability distribution of the 
internal energy (not shown). 

A similar dimensionality reduction has also been found 
in the ground state of the fee AF Ising model. 12 ) In addi- 
tion, similar magnetization plateaus have also been found 
in both the fee 13 )' 14 ) and pyrochlore models. 2 ) In these 
systems, however, the plateau states appear above some 
finite critical value of the applied magnetic field. A pe- 
culiar feature of our inverse perovskite model is that the 
critical field is zero and the plateau phases for the pos- 
itive and negative values of h touch with each other at 
h = with the complete suppression of T c , as shown in 
Fig. 3(c). 

The macroscopic degeneracy in the case of J = J± 
and h = is also lifted by introducing an anisotropy 
in the interaction, namely, J ^ J±. First, let us con- 
sider the case of J > Jj_. In the limit of J±/J — > 0, 
the system is decomposed into independent 2D square 
lattices and disconnected sites. Therefore, the system 
undergoes the same phase transition as that in the 2D 
Ising model at T c /J = 2/sinh _1 l ~ 2.269 15 ) and the 
ordered phase is the PDAF state in which the AF order- 
ings occur independently in the xy square lattices. This 
X7/-PDAF state is stabilized at low T in the entire region 
of J > J± . The order parameter is the staggered magne- 
tization of each xy plane; m xy = 

where the summation is taken over the spins within one 
xy layer. The transition is of second order, and the crit- 
ical temperature T c is determined by the Binder param- 
eter 16 ) g xy — 1 — (^^ 2/ )/3(m 2 ?/ ) 2 . The results are plot- 
ted in the left half of Fig. 4. T c decreases as J±/J in- 
creases and goes continuously to zero when J±/J — > 
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Fig. 3. (a) Temperature dependence of the specific heat at /i/ J = 
0.5 (left) and 4.0 (right) for J = J ± . The data at /i = are 
shown for comparison, (b) Magnetization process at several tem- 
peratures, (c) Phase diagram at J = Jj_. PDAF denotes the 
partially-disordered AF phase. The lines are guides for the eyes. 
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Fig. 4. Phase diagram for the anisotropy in the interaction 
(J, «/_!_) at h = 0. The lines are guides for the eyes. Schematic 
pictures of xy-PDAF state (left) and 3D Ferri state (right) are 
also shown. Open (filled) arrows denote up (down) spins, and 
balls represent disordered spins. 



1. We also examined the universality class of this xy- 
PDAF transition by the finite-size scaling (FSS) analy- 
sis; (m%. y ) = L~ 2 $l v f(tL x l v ), where / is a scaling func- 
tion and t = T/T c — 1 is the reduced temperature. We 
confirmed that the transition belongs to the 2D Ising 
universality with the exponents {3 = 1/8 and v — l. 17 ) 

On the other hand, on the opposite side, i.e., J < Jj_, 
we find a phase transition to the ferrimagnetic (Ferri) 
3D-LRO state, as shown schematically in the right in- 
set of Fig. 4. The order parameter is defined by mj = 
5^0^,1 + $v,2 — Sv,3)/3L 3 . The Ferri ground state 
has twofold degeneracy, rrtf = ±1/3 corresponding to 
the spin configurations (S^,i, 5^,2, ^,3) = (1, 1, — 1) and 
(—1,-1,1), respectively. In this case also, the transi- 
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Fig. 5. (a) Phase diagram of the inverse perovskite AF Ising 
model in the T-h-(J, J±) space, (b) h-T phase diagrams for sev- 
eral values of J±/J for J > J±. (c) (J, J±)-T phase diagram at 
h = 4max(J, Jj_). The lines are guides for the eyes, and TCP 
denotes the tricritical point. 



tion is of the second order, and the values of T c ob- 
tained by the Binder analysis are plotted in the right 
half of Fig. 4. T c decreases as J/ J± increases and goes 
continuously to zero when J/J± — > 1, as in the case 
of J > J±. We also confirmed that the transition be- 
longs to the 3D Ising universality class by the FSS anal- 
ysis. For instance, in the limit of J/J± — > 0, we ob- 
tained T c /J± = 4.0828(1) and the critical exponents 
/3 = 0.32(1) and v = 0.625(7), which are consistent 
with the 3D Ising exponents, (3 = 0.325 and v = 0.63, 18 ) 
within the statistical errors. 

The effect of anisotropy J j J± at h = is summarized 
as follows: The anisotropy lifts the macroscopic degen- 
eracy and induces the phase transition at finite T; the 
xy-PDAF state appears for J > J± and the 3D Ferri 
state for J < Jj_ . The two different ordered phases touch 
with each other at the isotropic point J = J±. Here, T c 
goes continuously to zero on both the sides. The univer- 
sality class of the transition is different between the cases 
of J > J_l and J < J±; the 2D Ising one in the former 
due to the dimensionality reduction and the 3D Ising one 
in the latter. These properties are contrasted to those of 
the fee system. In the isotropic fee model, a first-order 
phase transition takes place at finite T, namely, there is 
no suppression ofT c . 12 )' 19 ) It was also proposed that the 
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critical exponents change continuously 20 ) when a similar 
anisotropy in the interactions is introduced. 

On the basis of the above analyses obtained from 
applying magnetic field h or from introducing the 
anisotropy J ^ J± , the peculiarity of the isotropic model 
with J = J_l and h = is revealed: The model has a 
critical point at T — in the sense that the transition 
temperatures are suppressed down continuously to zero 
in both the cases of \h\ — > and J/ J± — > 1. 

To further clarify the behavior around the critical 
point, we investigate the model by changing both the 
h and J/ J±_ values. The phase diagram in a wide range 
of parameters is shown in Fig. 5(a). When we decrease 
J±/J from 1, two PDAF phases in the case of J = J± in 
Fig. 3(c) gradually merge into one, as shown in Fig. 5(b). 
On the other hand, in the region of J < Jj_, the magnetic 
field is a symmetry-breaking field of the Ferri ordering 
at h = 0, and therefore, the finite-T phase transition 
disappears at ft ^ 0. Hence, at finite ft, the PDAF or- 
dered phase suddenly terminates at J = J± . This behav- 
ior is studied in detail by changing J/J±. For instance, 
Fig. 5(c) shows the result at ft = 4max(J, Jj_). Since the 
phase transition at the isotropic case J = J± is of the 
first order as discussed earlier and since no transition 
exists at finite T for J < J_l, we have a discontinuous 
phase boundary at J = J± below T c . This discontinuous 
boundary slightly extends to J > J± and terminates at 
a tricritical point beyond which the transition becomes 
continuous, as shown in Fig. 5(c). The tricritical behav- 
ior appears for < \h\/J < 8, resulting in the tricritical 
lines which start from (Jj_/ J, ft/ J,T/ J) = (1,0,0) and 
end at (1,±8,0), as shown in Fig. 5(a). Note that the 
tricritical lines are not in the J = J± plane, but in a 
very narrow region, 0.95 < J±/ J < 1. 

All the features above are shown in Fig. 5(a) with the 
ground-state phase diagram which is easily obtained by 
minimizing the energy of one octahedron. The figure in- 
cludes the PDAF phase for J > J± (blue region), the dis- 
continuous surface in the ft = plane for J < J± (purple 
plane), and the tricritical lines (red curves) surrounding 
the discontinuous phase boundaries at the edges of the 
PDAF state (orange surfaces). This global phase diagram 
clearly demonstrates the peculiarity of the critical point 
at J = J_l and ft = 0, where all the critical temper- 
atures, for both first- and second-order transitions, are 
suppressed down continuously to zero. 

The absence of LRO at the critical point is understood 
as follows. When we start from a PDAF ordered state, 
we can generate an anti-phase domain in the AF-ordered 
square lattice without any energy cost because of the 
partial disorder. This leads to the destruction of PDAF 
ordering and to the formation of a degenerate ground- 
state manifold. The situation is similar when we start 
from the Ferri 3D-LRO state. These anti-phase domains 
cost a finite energy when we introduce the magnetic field 
ft and/or anisotropy J/J_\_ and hence finite-T phase tran- 
sitions immediately emerge around the critical point. 

In summary, by using the exchange Monte Carlo 
method, the phase diagram of the AF Ising model on 
the inverse perovskite lattice has been clarified in the 
parameter space of the uniform magnetic field and the 



anisotropy in the interaction. It has revealed an unusual 
character of the T = critical point, which drags the 
surrounding critical lines and surfaces at T ^ toward 
absolute zero temperature. Neither fee nor pyrochlore AF 
Ising model exhibits such a behavior. 

In the fee system at the isotropic point, there is no sup- 
pression of T c and a first-order transition to the 2D-LRO 
state takes place. 12 )' 19 ) In the isotropic pyrochlore sys- 
tem, although the ground state has macroscopic degener- 
acy as that in the present system, the degeneracy remains 
until a finite magnetic field is applied (|ft| < 2 J). 2 ) 

The pyrochlore system has attracted interest due to 
the degenerate ground state manifold, and the effects 
of the Heisenberg-type interaction and quantum fluc- 
tuations have been examined intensively in the search 
for an "order from disorder" phenomenon or a spin- 
liquid state. Since the degenerate ground state in the in- 
verse perovskite model is at a more unusual critical state 
with competing orders, it would be more interesting and 
worthwhile to study such effects in the inverse perovskite 
system. In particular, quantum fluctuations may easily 
destroy the orderings barely stabilized around the T = 
classical critical point, which may open possibilities for 
quantum spin-liquid phases or exotic quantum orders. 
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